Conductance fluctuations and field asymmetry of rectification in graphene 
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We investigate conductance fluctuations as a function of carrier density n and magnetic field in 
diffusive mesoscopic samples made from monolayer and bilayer graphene. We show that the fluctua- 
tions' correlation energy and field, which are functions of the diffusion coefficient, have fundamentally 
different variations with n, illustrating the contrast between massive and massless carriers. The field 
dependent fluctuations are nearly independent of n, but the n-dependent fluctuations are not uni- 
versal and are largest at the charge neutrality point. We also measure the second order conductance 
fluctuations (mesoscopic rectification). Its field asymmetry, due to electron-electron interaction, 
decays with conductance, as predicted for diffusive systems. 

PACS numbers: 



Reproducible conductance fluctuations (CF) are one of 
the most striking signature of phase coherent transport 
[1] . The conductance of a mesoscopic sample results from 
interference between all wave packets traversing the sam- 
ple. This interference pattern is sensitive to variations in 
disorder configuration, Fermi energy or magnetic flux, 
leading to reproducible CF as one of these parameters is 
changed. In diffusive or chaotic systems the CF ampli- 
tude has been shown to be universal [IHS] and ergodic, 
i.e. independent of the mechanism of phase randomiza- 
tion (magnetic field, Fermi energy, configuration of im- 
purities for diffusive systems, sample shape for ballistic 
systems). The CF amplitude is of the order of e^//i, with 
a coefficient which only depends on the symmetry class of 
the mesoscopic system. The typical correlation energies 
and fields B^^ of the fluctuations depend upon the 
typical time Tint and area Aint over which interference 
occur: = filnnt and = (f>o Mint, with $o = 
[2] . CF have been extensively investigated in metallic and 
semiconducting systems [U [3] . The recently discovered 
graphene [1] provides a unique system in which the Fermi 
energy and diffusion constant can be tuned at will, over a 
broad carrier density range extending from hole to elec- 
tron metallic conduction. Theoretical simulations of CF 
in graphene suggest a possible enhancement of the fluc- 
tuation amplitude with respect to standard mesoscopic 
samples, depending on the strength or nature of disorder 
(intervalley scattering) [SHZj. On the experimental side 
CF have been reported by several groups [8i413|. But 
to our knowledge the present work is the flrst complete 
investigation of their correlations and amplitudes as a 
function of Fermi energy and magnetic field, for both 
monolayer (ML) and bilayer (BL) graphene. The im- 
portance of the comparison lies in the fact that whereas 
ML and BL have similar resistivities and thus mean free 
paths (see Fig. [s]), the (massless and massive) carriers 
have different velocities because of the different disper- 
sion relations in these two materials. Thus the diffusion 
constants and therefore correlation energies and fields 
will have different carrier density dependences, provid- 
ing a powerful test of the applicability of the theory of 



mesoscopic fluctuations in those systems. We flnd that 
the variations with carrier density n of the correlation 
field and energy are well related to those of the diffusion 
coefficient. We also find, in contrast with [TO], that the 
amplitude of n-dependent fluctuations are largest near 
the charge neutrality point. 

We also measure the second order conductance fluctu- 
ations j21| which, unlike linear conductance, are a probe 
of Coulomb interaction and screening, an important issue 
in graphene. Second order CF are inherent to systems 
lacking spatial inversion symmetry (because of random 
disorder in diffusive systems or geometry in ballistic sys- 
tems). They stem from current-induced changes in the 
carrier density, which in turn, via Coulomb interaction, 
modify the electrostatic potential landscape, thereby in- 
ducing a current-dependent, or second order, change in 
the conductance of mesoscopic samples. Unlike the flrst 
order conductance which is even in field, the second order 
conductance has an odd part in field which was calculated 
for ballistic [TB] and diffusive [T7] systems. Whereas this 
mesoscopic rectification was experimentally investigated 
in ballistic mesoscopic systems (GaAs/GaAlAs quantum 
dots, Aharonov Bohm rings, carbon nanotubes HHHODj 
we provide in this letter the first measurement of second 
order CF in a diffusive system. We find, in qualitative 
agreement with theoretical predictions, that the odd in 
field rectification decreases with conductance. 

Both ML and BL graphene samples (the nature of 
which was confirmed by Raman spectroscopy) were exfo- 
liated and deposited onto doped silicon substrates with 
a 285 nm thick oxide. The electrodes were fabricated 
by electron beam lithography and sputter deposition of 
40 nm thick palladium, providing low contact resistances 
(200 n for the ML and 20 Vt for the BL). The dimen- 
sions are W — 2.7 /im, L — 0.8 jim and W = 4.8 /^to, 
L = 0.7 jim for the ML and BL, respectively (Fig. |3|. 
The carrier density n is controlled for both systems by 
the voltage Vg applied to the doped silicon back gate, via 
n ~ kp/n — C{Vg — Vq), with C the gate capacitance 
per unit area, kp the Fermi wave vector, and Vb the gate 
voltage of the charge neutrality point. The transport pa- 
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rameters of both samples were determined via classical 
magnetoresistance measurement, described in [22]. We 
found that the transport mean free path Itr of both ML 
and BL varies between 25 and 100 nm over the 60 V gate- 
voltage range probed. Thus Itr is much smaller than the 
sample dimensions, so that both samples are in the diffu- 
sive regime. Two terminal resistance measurements were 
performed at 60 mK in a dilution refrigerator with resis- 
tive lines filtered at room temperature. The ac current 
amplitude Iq was adjusted between 5 and 50 nA, de- 
pending on the sample resistance, to avoid heating while 
optimizing the signal to noise ratio. The first and second 
harmonics Vi and V2 of the ac voltage were measured 
with a low noise voltage preamplifier and a lock-in am- 
plifier. The first and second order conductances were de- 
termined via Gi = lo/Vi and G2 = 2^2 /o/V? [H]- The 
magnetic field was limited to less than 0.4 T, so that the 
contribution of Shubnikov de Haas oscillations is negli- 
gible. The CF of Gi were measured as a function of B 
and Vg, for both the ML and BL (see Fig. [T]). The Vg 
dependent fluctuations were analyzed in 3 V-wide win- 
dows, over which the density n can be considered practi- 
cally constant. The fluctuations of Gi, even in B as ex- 
pected in a two probe measurement, were characterized 
by their correlation functions and their amplitude SGi, 
defined as the square root of their variance. The his- 
togram of these fluctuations is Gaussian over the whole 
parameter range. The correlation gate voltage and fleld 
were determined in the following way. We calculated the 
Fourier power spectrum of the fluctuations as a function 
of k, the conjugate variable of B or Vg, for each set of 
fleld or gate voltage data corresponding to a given aver- 
age carrier density. This spectrum is the Fourier trans- 
form of the correlation function of the fluctuations. Each 
dataset yielded exponential-like functions which corre- 
spond to a lorentzian correlation function in the direct 
space (Fig. [2^). The decay exp{—kxc) at low k yields 
the correlation gate voltage or field Xc for each dataset. 
The correlation energy Ec is related to the Vg correlation 
scale Vc by = VchvF\/T^C ({/^eVg) for the ML and 
7i^7rGT4/(2eme//) for the BL, with TOe// = O.OSmg the 
(constant) low energy effective mass in the BL. These 
relations are incorrect near the charge neutrality point 
(CNP) because of density inhomogeneities, so that we 
deduced E^ from Vc via kp{Vg) determined in [^ . 

In Figj2] we compare the variations of Ec and Be 
to the predictions ^ E^ — TiD/L^^^ and B^ = 
^o/{LminWrmn), whcrc Lmm = min{L,LT,L^) and 
Wmin = min{W, Lt, L^) correspond to the typical longi- 
tudinal and transverse length of interfering trajectories; 
Lt = ^TiD /ksT and are respectively the thermal 
and the phase coherence length. We first determine the 
diffusion coefficient D = a / {e^v{Ep)) — vplfr/'i, where 
a is the conductivity and v{Ep) is the density of states, 
vp is constant (10^ m/s) for the ML and is hkp/meff 
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FIG. 1: CF as a function of magnetic field (left) and gate 
voltage (right). A linear fit was subtracted to remove the 
average conductance. Thick red (resp. thin blue) curves are 
measured around Vg — Vo, the charge neutrality point (CNP) 
(resp. at high carrier density, corresponding to Vg — Vb = 
40.5 V for the ML and Vg-Vo = 15V for the BL). The CNP 
is at Vo = 15.5 V for the ML and Vo = -6 V for the BL. It is 
clear (b,d) that the amplitude of Vg-dependent CF is larger 
at the CNP than at high density, both for the ML and the 
BL. Panel d also shows the reproducibility of the CF, via two 
nearly indistinguishable successive gate voltage sweeps (open 
dots and red line). 



for the BL. Both Itr and kp were extracted in [22]. We 
find that at 60 mK Lt varies between 0.7 and 1.8 /xm, so 
that L < Lt < W for both samples over the entire den- 
sity range investigated. Thus Wmin — Lt and Lmin — L. 
Consequently Ec is expected to vary like the Thouless en- 
ergy Etk = fi/TD^ hD/L"^, and Be like ^q/{LLt) [H]. 
As shown in Fig. 2, this is indeed what is found: the 
correlation energy deduced from the experimental data 
follows Ec = 5 Eth, and the correlation field follows 
Be ^ 6 ^o/[LLt) for the ML and = 4 ^q/{LLt) 
for the BL. In particular, we find the expected linear de- 
pendence of Ec with D, which itself varies like Vg for the 
BL and ^/Vg for the ML. This stems from the fact that in 
both samples the conductance varies linearly in n cx Vg 
within logarithmic corrections [22] . and that v{Ep) is in- 
dependent of Vg for the BL and varies like ^fVg for the 
ML. Similarly Be oc ^q/LLt is expected to vary like 
D^^/^ cx Vg'^^'^ and Vg~^^^, as observed respectively for 
the BL and ML. The numerical factors may be explained 
by the samples' aspect ratios, but more theory is needed 
to assess this point. 

We now discuss the CF amplitude. As is visible in fig- 
ures [T|d and [iji, the Vg-dependent CF are greater near 
the CNP. This is quantified by the fiuctuation amplitude 
6G1 plotted as a function of Vg in Fig. |3] In contrast, the 
-B-dependent fluctuations do not significantly depend on 
the density n. This non ergodicity of the fluctuations in 
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FIG. 2: Correlation energies and fields of CF. (a) Fourier 
power spectrum of CF at different Vg in a semi-log scale. T4 
is extracted from the exponential fit at small AVg~^. (b) 
and (d): EdVg) for the ML and BL. Continuous lines are 
5ETh, with D determined from previous magnetoresistance 
measurements (see text), (c): Bc{Vg). Continuous lines are 
= 6<1>o/(-LtI') for the ML and 4$o/(iTi) for the BL, 
with D similarly determined. Vertical error bars represent 
uncertainty in fit parameters, and horizontal bars the 3 V 
range used to calculate correlation energies. 



graphene may be a consequence of the spatial inhomo- 
geneities of n close to the CNP. Indeed, in a good con- 
ductor (large g) , changing the Fermi energy is equivalent 
to changing the disorder configuration, and induces CF 
of order /h. In graphene on the other hand, it has been 
shown [321 that near the CNP the system breaks into con- 
ducting puddles of electrons and holes, and that trans- 
port takes place along an intricate percolating network 
of these n and p-type regions iGj . Thus, as pointed out in 
[5], the change in configuration induced by a change of 
Vg in this region may induce a much larger variation of 
conductance than a change in the disorder configuration 
in a good conductor. In contrast, the magnetic field does 
not affect the network but only the phases of the wave 
functions, which explains why the fluctuations with B 
are independent of gate voltage jS]. Note that these re- 
sults are in contradiction with those of [TUI who found 
a smaller CF amplitude close to the CNP. This discrep- 
ancy could be due to a different nature of disorder in the 
samples ^ . For the ML, we find the i3-dependent CF 
amplitude to be 0.4 ± 0.1 e^/h, independent of Vg. The- 
ory [3j predicts 6G = O.Ty^W/L e^/h when the distance 
between electrodes L is smaller than Lt- This yields 
SG = 1.2 e^/h for the ML, which is three times the value 
measured in this work (Fig. [3|. 

We now turn to second-order non-linear CF in mono- 
layer graphene. As mentioned in the introduction, cur- 
rent through a mesoscopic sample induces charge ac- 
cumulation around impurities and sample edges. This 
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FIG. 3: Comparison between resistance and CF. The Vg- 
dependent fluctuations (filled symbols) are larger near the 
CNP for both the ML and BL (panels c and d). Error bars 
are the standard deviation of the difference of two differ- 
ent sweeps. Dashed lines are guides for the eyes. The Vg- 
dependent fluctuation amplitude changes with Vg in a quali- 
tatively similar way as the resistance. The _B-dependent CF 
amplitude (open triangles) doesn't change much with Vg. 



bias induced change of electronic density Jn(r, V) mod- 
ifies (via Coulomb interactions) the electrostatic poten- 
tial landscape throughout the sample by SUdis {lint , V) , 
where quantifies e-e interactions jl5j . As a result, 
CF acquire a non linear bias dependent contribution : 
G = G[Ueq]+5G[5Udis{nint, V)], where Ueq is the poten- 
tial in zero bias. Gi — G\Ueq{r) is the linear conduc- 
tance, and G2 — {dG[5Udis^V,iint\/dV)v=o is the sec- 
ond order conductance. In a magnetic field -B, 5n{r, V) 
and SUdisilint, V) contain a component which is even, as 
well as another which is odd in B, which can be viewed 
as a mesoscopic Hall voltage. Contrary to Gi which is 
even in i3, in the presence of interactions G2 has a com- 
ponent which is odd in G^'^. The ratio between the 
even (symmetric) variance <5Gf and the odd (antisym- 
metric) variance JG^'^ of G2 is predicted to be indepen- 
dent of conductance in ballistic systems, but in diffusive 



systems it should vary like -j^s~ 
ductance in units of /h [T71 1 



f^, with g the con- 



Figure [4] shows Gf 
and G^'^ at two different carrier densities (close and far 
from the CNP), as well as the ratio of the odd to even 

SqAS 

amplitude r = . We find that r significantly de- 

creases as carrier density increases, in contrast to bal- 
listic GaAs/GaAlAs rings [20], where r was nearly inde- 
pendent of the dimensionless conductance g. To compare 
to the predictions for diffusive systems, r = lint/ 9, we 



plot in Fig. 
cording to 



r(Vg), and find that r « (7 ± Ac- 



£], the interaction constant 7i„t is related 
to the geometrical capacitance C and electro chemical 
capacitance = (l/(j^e^) + I/G)^^ by "fi^t — Cfj_/G. 
Thus lint = 1/(1 + Cjive^)) [TB]. Perfect screening (in 
materials with large densities of states) corresponds to 
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FIG. 4: Second order conductance fluctuations for tlie ML, 
near (a) and far (b) from the CNP. Tfie antisymmetric part 
{G2{H) - G2{-H))/2 is largest at the CNP. (c)The 



ratio r — 56*2^* /^Cf is compared to 7/g. Error bars are de- 
duced from the difference of G2{H) between two consecutive 
scans, which is mainly due to telegraphic noise in the sample. 



7mt = 1, and no screening to 7i„i = 0. In graphene 
7mt ~ 1 over the entire n range investigated, within less 
than 4/1000 : screening is strong, even close to the CNP 
where jmt = 0.996. The factor 7, compared to 1 expected 
for L — W, may be due to the large aspect ratio W/L of 
the sample, which is known to increase SGi entering in 
the calculation of r ^ , but there are no calculations for 
our geometry. 

In conclusion we have shown that mesoscopic graphene 
samples exhibit conductance fluctuations which Fermi 
energy- and i3— dependent correlation functions can be 
described by theoretical predictions for diffusive systems 
over a wide range of carrier concentration, for both mono- 
layer and bilayer. The different behaviors of the correla- 
tion energy and fields are intimately related to the funda- 
mentally different dispersion relations of both systems. A 
significant increase of the amplitude of the Fermi energy- 
dependent fluctuations is observed close to the neutrality 
point, whereas the i?-dependent fluctuation amplitude 
is nearly constant over the entire carrier density range. 
This non-ergodicity of fluctuations may be attributed to 
the particular disorder due to electron and hole puddles 
in graphene near the Charge Neutrality Point. Finally, 
we have measured the second order non-linear conduc- 
tance. We have exploited the tunability of graphene's 



conductance to find that its field asymmetry decreases 
with g, in agreement with theoretical predictions for dif- 
fusive systems. This indicates strongly screened electron- 
electron interactions in graphene. 
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